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Abstract. We study the unoccupied electronic structure of the spin-i quantum 
magnet TiOCl by x-ray absorption near-edge spectroscopy (XANES) at the Ti L 
and O K edges. Data is acquired both in total electron and fluorescence yield mode 
(TEY and FY, respectively). While only the latter allows to access the unconventional 
low-temperature spin-Peierls (SP) phase of TiOCl the signal is found to suffer from 
significant self-absorption in this case. Nevertheless, we conclude from FY data that 
effects of the SP distortion on the electronic structure are absent within experimental 
accuracy. The similarity of room-temperature FY and TEY data, the latter not 
being obscured by self-absorption, allows us to use TEY spectra for comparison with 
simulations. These are performed by cluster calculations in D4/J and D2/1 symmetries 
using two different codes. We extract values of the crystal-field splitting (CFS) and 
parametrize our results in often found notations by Slater, Racah and Butler. In all 
cases, good agreement with published values from other studies is found. 
(Some figures in this article are in colour only in the electronic version) 
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1. Introduction 

Low-dimensional quantum spin systems often have a rich phase diagram, especially when 
they involve magnetic frustration. If strong electronic correlations are also important 
due to partially filled d or f shells there will be strong interplay between spin, charge, 
orbital, and lattice degrees of freedom. In such systems a method like x-ray absorption 
spectroscopy (XAS) although it measures only the unoccupied density of states (DOS) 
thus gives valuable information about all degrees of freedom, not only the electronic ones. 
Figure [1] shows the crystal structure of the layered transition metal oxyhalide TiOCl. It 
has a spin-Peierls (SP) ground state, i.e. a dimerized phase due to dominant coupling 
between lattice and spin [1] . This state is reached, however, in an unconventional fashion 
by two consecutive phase transitions [21 |3l H] . That is, upon lowering the temperature a 
kink in the magnetic susceptibility is observed at Tc2 = 91 K, and incommensurate peaks 
are found in x-ray diffraction (XRD). At = 67K the susceptibility has a sharp drop 
and shows hysteresis, concomitant with a lock-in of the XRD peaks at commensurate 
(hkl) positions, i.e. the dimerization of Ti ions is completed [5]. This succession of a 
second-order and first-order phase transition signals the non-canonical character of the 
SP transition in TiOCl. While the electronic structure at elevated temperatures (room 
temperature and above) has been investigated intensively both by theory O O El |9l [10] 
and experiment [HI |12l |T3l [HI [15] , the low-temperature phases are not accessible by 
photoemission spectroscopy (PES). TiOCl is a 3d^ Mott insulator, and below room 
temperature its strongly insulating character leads to excessive charging and a loss of 
all spectral features in PES. Fluorescence yield (FY) x-ray absorption spectroscopy can 
overcome this problem since the detected photons are not susceptible to sample charging. 
As an analysis of FY spectra introduces further complications like e.g. self-absorption 
we have used total electron yield (TEY) data more extensively, as will be discussed 
below with appropriate caveats. 

The potential of XAS at the oxygen K and transition metal (TM) L edges to 
investigate the electronic structure of strongly correlated electron systems has been 
shown for many materials in the past two decades [161 [13 [HI [IS [201 [HI [221 123] . A simple 
model case is when the TM ion has a 2p^3d^ valence configuration and is octahedrally 
coordinated by six oxygen ligands (cf. Fig. [1]). This gives rise to crystal- field (CF) 
splitting of the 3d orbitals into triply degenerate t2g and doubly degenerate Cg levels 
due to the local O^ symmetry. Hybridization of the TM 3d shell with the O 2p ligand 
orbitals allows an investigation of the unoccupied DOS of the TM ions also from the O 
K edge spectra. An advantage of the O K edge spectrum is the absence of multiplet 
splitting which would lead to rather complicated shapes of the absorption spectra, as 
seen at the TM L edge. The latter, however, can be tackled by local calculations on single 
clusters (i.e. (TM)06) using the configuration interaction cluster model with full atomic 
multiplet theory. This is highly powerful when simulating L edge spectra because one 
can include the initial and final states as well as spin-orbit coupling, multiplet splitting, 
and charge transfer effects on an equal footing. 
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Figure 1. Crystal structure of TiOCl. Left panel: Two double layers separated 
by van-der-Waals gaps. Center panel: Comparison of an ideal octahedral coordination 
(O/i symmetry; top) and the actual orthorliombic coordination of Ti ions in TiOCl {C2v 
symmetry; bottom). Right panel; Projected view along the c axis of the Ti sublattice. 
Top: Undimerized high-T structure. Bottom: Dimerized spin-Peierls phase. Thick 
blue bonds indicate = singlet pairs. 

The Ti ions in TiOCl lie in a strongly distorted O4CI2 octahedron (cf. Fig. [1]), 
resulting in five well separated CF levels and a quenching of the orbital degree of freedom 
[HI [lU [21]. The local point-group at the Ti site is C2V, with the C2 axis parallel to the 
crystal c direction, pointing in between two CI ligands. In this symmetry, the dx2_y2 and 
the orbital belong to the same irreducible representation (ai) and therefore can mix. 
It is found that the CF eigenstates in TiOCl are close to a (i^.2_j^2-like orbital pointing 
in the c and b direction for the ground state and a (i^2-like orbital pointing in the a 
direction for a highly excited state within the commonly used reference frame (x = b, 
y = c, z = a) . Accordingly, the eigenstates can be labeled (in ascending energetic order) 
ai (ife2_c2 = dx2-y'2, a2 dab = d^z, bi dac = dyz, 62 c?bc = dxy, and ai da^ = dz^, i.e. the only 
electron in the 3d shell has dx2-y2 character. This provides a preferred one-dimensional 
hopping path along b due to direct overlap of electron clouds. As the XAS final state 
now has two 3d electrons, multiplet mixing is refiected both mO K edge and Ti L edge 
spectra. 
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2. Experimental and technical details 

Single crystals of TiOCl were prepared by chemical vapor transport from TiCls and 
Ti02 [2S]- Typical sample dimensions were 1 x 3 x 0.2 mm'^. They were cleaved in situ 
at pressures below 1 ■ 10~^mbar using Scotch tape to expose fresh surfaces. Experiments 
were performed with the MUSTANG endstation at beamline PM3 at BESSY (Berlin, 
Germany) with parallel recording of the FY and TEY signals. For the latter a Keithley 
Nanoampere meter was used to measure the drain current from the sample which was 
mounted in an isolated fashion. Absorption features (in particular carbon derived) 
of the plane grating monochromator were accounted for by dividing the TEY and 
FY signals by a simultaneously recorded spectrum from a gold mesh. The overall 
resolution amounted to 100 meV full width at half maximum (FWHM) and the base 
pressure in the experimental chamber was below 5 ■ 10~^*^ mbar. Data was recorded in 
normal-incidence geometry, and the effective polarization was changed by rotating the 
sample stage around the beam axis. The manipulator allowed for cooling the sample 
down to at least 80 K (checked by a thermocouple on the sample holder), i.e. well 
below the onset of the incommensurate lattice distortion at Tc2 = 91 K. We compare 
different codes for our cluster calculations: crystal-field theory (CFT) calculations were 
performed using program codes originally developed by R. D. Cowan [26j, P. H. Butler 
[27], and B. T. Thole and F. M. F. de Groot [H]. As the desired symmetry is reached 
by successive symmetry reduction via different point groups it becomes increasingly 
difficult to define correct (and meaningful) CF parameters and operators for given 
polarizations. Thus, comparatively high tetrahedral symmetry {D^h) with a distortion 
along the crystallographic a axis was used to simulate XAS spectra of TiOCl. Ligand- 
field multiplet (LFM) calculations in D2/1 symmetry were performed with a code by M. 
W. Haverkort. The ligand-field theory parameters are calculated using Wannier Orbitals 
obtained from a density-functional theory (DFT) calculation done in the local-density 
approximation (LDA) with the Stuttgart LMTO code |28j . 

3. Results and discussion 

Figure [2] shows a complete set of O edge FY and TEY spectra taken at room 
temperature and at 80 K for linear polarization both along the a and h axis. From 
the FY data one sees that no significant changes in the electronic structure take 
place between the high-temperature phase and the incommensurate intermediate phase 
(Tci < 80 K < Tc2). Since this observation is resolution-limited, changes in the electronic 
structure occur obviously only on the energy scale of some 10 meV or less. This supports 
the assumption that the charge degrees of freedom are barely involved in the spin- 
Peierls transition in this oxyhalide. From the very good correspondence between room- 
temperature FY and TEY spectra [cf. panel (b)] we conclude that, despite a reduced 
probing depth, TEY spectra reflect bulk properties of TiOCl and can be used for a 
detailed analysis. Thus, one circumvents possible problems from self- absorption which 




Figure 2. O K edge absorption spectra of TiOCl. Spectra with polarization along b 
are vertically offset for clarity, (a) FY spectra at room temperature (yellow diamonds) 
and 80 K (blue squares). The ranges of predominantly Ti 3d [hv < 537 eV) and Ti 
4s/p [hv > 537 eV) bands are indicated at the bottom. For 3d derived structures, the 
peaks of mainly t2g and mixed t2g and Cg character are marked (see text for details) . 
(b) Comparison between FY (yellow diamonds) and TEY (red squares) spectra at 
T = 300 K. FY spectra obviously suffer from self-absorption; arrows mark prominent 
examples. Note that from the room-temperature TEY data along a, pre-edge data 
points between 530.0 and 531. 7eV have been excluded as bad data points. 



appear e.g. in the FY spectrum with polarization along b, where the intensity of the peak 
at lowest energy appears suppressed. The observed absorption features can roughly be 
assigned to O 2p states mixed with Ti 3d {Eph ~ 531 — 537 eV) or 4s/4p {Eph > 537 eV) 
orbitals, respectively, and the double-peak structure in the low-energy part comes mostly 
from the (quasi-) octahedral coordination of the TM ion. In an ideal Oh symmetry this 
peak separation would directly reflect the t2g-eg CF splitting lODq [18]. Although 
TiOCl belongs to the point group with considerably lower symmetry, the broadness 
of features as well as the illustrative value of the terms from 0^ symmetry justify usage 
of this nomenclature for qualitative arguments. 

The final state of the XAS process in TiOCl involves two electrons in the d shell, 
which basically allows for triplet and singlet configurations with both electrons in t2g 
orbitals (t^^; ^Ti for S = 1; ^T2, ^E, and ^Ai for S = 0), one in a t2g and one in an 
eg orbital [tlge]^; triplet states are ^T2 and ^Ti), as well as t^g^g configurations. The 
latter as well as t^gCg singlet configurations can be neglected, since they lie in energy 
well above all others and thus certainly do not contribute to the first two peaks. Table [H 
gives the basic energetics of these configuration^ in terms of U, J and lODq [29], as well 
as values calculated by Macovez et al. who applied an embedded cluster approach to O 

t ^2g'2g singlet configurations arc included only for completeness. 



Unoccupied electronic structure of TiOCl 



6 



Table 1. Possible configurations, term symbols and energy eigenvalues of a system 
in Oh symmetry {t2g threefold degenerate). The definitions of U and J in terms of 
Racah parameters can be found in table lAll In the second-to-last column, energy 
ranges (in eV) relative to the ground state {^Ti) from an embedded cluster approach 
to O K edge spectra of TiOBr are given (taken from [SO])- The last column gives 
energies (in eV) relative to the ground state, calculated using the best-fit parameters 
to the Ti L edge spectra from CFT calculations in D/^h symmetry (see text for details) . 



configuration 


term(s) 




eigenenergy'^ 


range in TiOBr 


best-fit energy 


t\g triplet 


^Ti 




U-3J3 







t\g singlet 










1.44 


tig singlet 


'A, 




U + 2J3 


0.78-1.28 


3.40 


*2geg triplet 




U - 


- 3 J4 + IQDq 


1.32-1.61 


1.61 - 2.52 


*2geg triplet 




U - 


- 3 Ji + lODq 




1.90 - 2.82 


tlgf^l singlet 




u 


- Ji + lODq 




1.65 - 3.25 


*2geg singlet 




u 


- Ji + lODq 




3.34 - 3.65 



^ Notation adopted from |29j . 

^ The absolute value using U and J is 3.04 eV. 



K edge spectra of TiOBr [30]. The energetic order found by these authors [30] is such 
that the triplet configuration ^Ti is at lowest energy, followed by the singlet ^Ai , and 
finally the and ^Ti t\ge]j triplet configurations. Due to symmetrically inequivalent 
O sites in the applied supercell, they could give only an energy range for the different 
configurations from their calculation, depending on the position of the induced core hole 
on a specific O site. It is well justified to assume that these results and assignments 
also apply to TiOCl. We thus identify the peak at Ep^ ~ 532 eV as the ^ig^g triplet 
configuration, while the one centered at Eph ~ 534 eV is a superposition of ^Ai singlet 
and triplet configurations. Calculating the eigenenergies from the relations given 
in the third column of table [T] either with U = 4.5 eV, J = = 0.7 eV from GGA [5] 
OT U = 5.2 eV, J = 0.72 eV from best fits within our CFT calculations (last column in 
table [H see below for details) shows how important the attractive force between core 
hole and excited electrons is in x-ray absorption: in neither case can the results from 
the embedded cluster approach (which takes this Coulomb interaction correctly into 
account) be reproduced by the simple d'^ energy scheme. Nevertheless, it is useful in 
situations where no core hole is present, as is evidenced by comparison to PES studies 
on electron-doped TiOCl [21] • 

Upon switching the light polarization between a and b, two major changes are 
observed. Firstly, the s/p hybridized bands change their shape significantly. However, 
due to the large number of bands involved here, a simple interpretation of this behaviour 
is not possible, but these features serve well as a cross-check of the sample orientation. 
Secondly, the intensity ratio of the two 3d peaks is changed, making the first peak at 
Eph ~ 532 eV higher than the one at Eph ^ 534 eV for polarization along b. However, 
one cannot use simple dipole matrix arguments (i.e. selection rules) to explain this 
behaviour due to the mixed multiplet character of each feature. Instead, a calculation 
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Figure 3. Ti L edge absorption spectra of TiOCl at T = 300 K with polarization 
(XAS) or momentum transfer (EELS) along the (a) a and (b) b axis. EELS data 
(blue triangles) provides the correct energy scale. Also shown are TEY measurements 
(yellow squares) and results from CFT (blue dotted curves) and LFM (red full curves) 
calculations in D4/J and T)2h symmetry, respectively (for details see text). The insets 
show corresponding FY spectra at room temperature (full yellow diamonds) and 80 K 
(open blue squares). Self- absorption effects appear even stronger than in O K edge 
spectra, cf. structures marked by arrows. 



of the Ti partial DOS e.g. using density-functional theory would be required. 

Figure [3] presents measurements of the Ti L edge for polarization along a and h 
from different methods. The EELS data was recorded with a dedicated spectrometer 
[32] using a primary beam energy of 172 keV and an energy resolution of AE ^ 270 meV 
at T = 300 K. The sample preparation for these measurements and further experimental 
details are described in [33]. The EELS energy scale was used to calibrate the photon 
energy scale for the XAS data. The insets show FY spectra both at T = 300 K and 
T = 80K, which appear to suffer even more from self-absorption than the O K edge 
spectra in figure [21 No differences can be observed between the high-T and low-T phases, 
similar to the O K edge. As discussed before, we used TEY data to compare with our 
simulations, as they show finer structures and need not be corrected for self-absorption. 
Results from calculations within CFT in D4/1 symmetry and from LFM theory in 
symmetry are also shown in the figure. All theoretical spectra have been broadened 
by a Gaussian with FWHM of 0.1 eV to account for the experimental resolution, and 
Lorentzians with FWHM of 0.4 eV and 0.2 eV for the L2 and L3 edges, respectively, to 
account for lifetime broadening. Different values for the two edges are justified by the 
fact that there is an additional Auger decay channel available at the L2 edge, leading to 
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Table 2. Best- fit input parameters for CFT caicuiations in symmetry to simulate 
Ti L edge TEY spectra, and resulting crystal-field, Racah and Hubbard parameters. 
The first four rows are the input parameters to the programs, and the last four rows 
have been calculated from these according to known relations, cf. table ETl [M] . The 
values given for the Slater integrals are the atomic parameters from Hartree-Fock, 
which were reduced to 80% of that value to give the best fit.'' Crystal-field splitting is 
given relative to the ground-state orbital. All values are given in eV. 



spin-orbit coupling 


2p 

3.8 ±0.4 


M 

0.0272 ±0.01 




direct Slater integrals 


^ pd 

5.581 


F2 
^ dd 

10.343 


F^ 
^ pd 

6.499 


exchange Slater integrals 


^ , 1 

3.991 


^pd 

2.268 




crystal-field integrals'^ 


X400 
4.9 


X420 

-0.05557 


X220 
1.29 


crystal-field parameters'' 


IQDq 
1.5 


Dt 

0.00343 


Ds 

-0.1543 


crystal-field splitting'^ 


0.48^ 


1.48 


6,^2 

2.08 


Racah parameters 


A 

3.559 


B 

0.101 


C 

0.413 


Hubbard parameters 


U 
5.6 


0.74 





^ It is well-known that this reduction mimics the Slater integrals in actual solids |23j . 
'^ Apply only to CFT calculations. 

^ These resemble the values observed in resonant inelastic x-ray scattering (RIXS) 
[24] . The latter were used as initial guesses for the crystal-field parameters in CFT 
calculations. In LFM calculations, they were fixed input instead of fit parameters. 
This is the weighted average of the dxz and dyz CF splitting seen in RIXS [24] . 

enhanced broadening [IB]. The best-fit parameters for the CFT calculations are given 
in table |2l Note that d^z and dyz are degenerate in 04/^ symmetry, and the value of the 
CF splitting given in the table thus can be understood as their weighted average. 

Comparing first the CFT calculations (blue dotted curves) with the experimental 
spectra one realizes that although the basic shape in terms of number of peaks, 
relative intensity and polarization dependence can be reproduced, they appear stretched 
energetically compared to experiment. Such stretching can be straightforwardly 
attributed to the omission of charge-transfer effects (i.e. hybridization with the ligands) 
when one looks at the spectra calculated within LFM (dotted curve). It is well 
known that these effects lead to a contraction of multiplet structures if the charge- 
transfer energy A (= 5.0 eV for the shown spectra) is positive, and produce additional 
small satellite structures [20]. Also, the degeneracy between d^z and dyz is lifted in 
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these calculations due to the lower symmetryl§| Nevertheless, transmuting the 
fit parameters from the different methods into the on-site Coulomb repulsion U and 
exchange interaction J (= J3 for the cases relevant to TiOCl [29]) one finds good 
agreement with typical values needed for reasonable results from density functional 
calculations [8]. 

In order to arrive at the LFM spectra, a tight-binding representation of the band 
structure on a basis of Ti-ci, 0-p and Cl-p Wannier orbitals was calculated using the 
A'^MTO approximation [36]. Starting from these Wannier orbitals we defined local 
correlated Ti orbitals and uncorrelated ligand orbitals. The ligand orbitals were obtained 
as a Ti-centered linear combination of the Wannier orbitals centered at the O and CI 
site. More details on this procedure can be found elsewhere [28]. For systems with a 
local cubic point group there is only a single ligand shell. In our case, the hybridization 
of the first ligand shell with the second turned out to be non-negligible and is included. 
For each d orbital there are two ligand orbitals, labeled L^^^ and L^'^\ The d shell has 
interaction with the L^^^ shell, which in turn interacts with the L*^^^ shell. The non- 
spherical part of the interaction found in this way is given by the following irreducible- 
representation-dependent one-particle Hamiltonian: 



H, 



0.2 



I 


dxz 


r(l) 
-l-'xz 


r(2) \ 
-I^xz ' 




dxz 


-0.461 


1.862 






r(l) 

J-/XZ 


1.862 


-0.808 


1.054 


\ 


r(2) 
-l-'xz 




1.054 


1.592 / 



(1) 



Hi 



b2 



/ 


dyz 


r{l) r{2) 
^yz ^yz 




dyz 


0.013 


2.060 




r(l) 
J-'yz 


2.060 


-0.399 1.084 








1.084 0.929 / 


/ 


dxy 


r(l) r(2) \ 
i-'xy i-'xy ' 




dxy 


0.325 


2.529 




r(l) 
J-'xy 


2.529 


1.075 1.509 




\ r^^) 




1.509 0.562 ) 





(2) 



(3) 



ai 





^X^—'ip- 


dz2 


^(1) 




r(2) 

^ ?/ 


4? 






-0.564 


-0.081 


1.528 


-0.123 










-0.081 


0.713 


-0.123 


3.273 








r(l) 


1.528 


-0.123 


-0.062 


-0.103 


1.083 


0.110 




z^ 


-0.123 


3.273 


-0.103 


0.194 


0.110 


0.910 




r(2) 

J-J n 






1.083 


0.110 


0.930 


-0.086 


v 


2^ 






0.110 


0.910 


-0.086 


0.908 



(4) 



§ In V)2h one can define two additional crystal- field parameters j3S|. While Du directly describes the 
splitting between dxz and dyz , the parameter Dv gives a mixing (and thus additional splitting) of (i^2 
and dxy orbitals. 
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The LMTO potential used to obtain the downfolded Wannier orbitals and tight- 
binding representation of the band structure was calculated using an atomic sphere 
approximation (ASA). Although reasonably accurate, the non-spherical part of the 
Madelung potential is missing within this potential. We included this as a set of ad-hoc 
parameters, fitted to the d — d excitation eigenenergies as obtained from the resonant 
inelastic x-ray spectra of TiOCl [21]. The resulting additional potentials are: 

dx^—y'^ dxz dxy d^^ 1 

-0.036 -0.156 0.144 -0.056 0.104 ^ ^ 

Within the LFM calculations the center of gravity of the d^ configuration is taken to be 
at the energy 0. The d'^L configuration has its average energy at A = 5.0 eV, whereby 
L can either be a hole in the first or the second ligand shell. The d^L^ configuration 
has its average energy at 2A + Udd with Udd = 4.0 eV, and so forth. 

4. Conclusions 

We study the unoccupied electronic structure of the frustrated quantum magnet TiOCl 
using Ti L and O K edge x-ray absorption near-edge spectroscopy. From total electron 
yield and fluorescence yield data we conclude that the low-temperature spin-Peierls 
distortion has no influence on the electronic structure within experimental accuracy. 
Comparing our data with cluster calculations in D4/1 and symmetries using different 
codes we can extract values for the crystal-fleld splitting which are in line with other 
work, and parametrize our results in widely found notations by Slater, Racah and Butler. 
Values of the extracted electronic parameters like U and J are in very good agreement 
with those from advanced methods like e.g. density-functional theory. 
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F4 



1 p2 

A 4 1 



^ — r 441 -r 

n _ 1 p2 5 p4 



F^ = Fo Fo A = Fo- 49F4 



normalized SCP F^ = 49F2 Fo B ^ F2 - 



2 ^2 D ~ r2 ~ or^ 

F* = 441^4 ^4 C = 35F4 



F" = A - |C Fo = A~lC A 

Racah F^ ^ 49^ _ Fa = S - iC B 

7Ti4 441 J? — i r' n 

b - 3g-G i^4 - 35O C 



u 








4F2 + 


36F4 


A + 4,B^ 


-3C 


U' = 


U - 2Jh 


49 441^ 




2F2 - 


4F4 


A- 


25 H 


-C 


Jl 




35 p4 
441^ 


35F4 






C 








Jh) 


A f2 , ^ p4 
49 ^ 441 


3F2- 


I-2OF4 




SB 


+ C 




■h 




A P2 , JJ. p4 
49-^ ^ 441^ 


4F2- 


I-I5F4 




AB 


+ C 
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6Dq + 2Ds - Dt 


;^^400 - ;^^420 - -^^220 




aig 


6Dq - 2Ds - 6Dt 


■^Xioo + :^X420 + ■^X220 


dxy 


b2g 


-ADq + 2Ds - Dt 


3730^^-"^" ' 3^42^420 ^70^220 




eg 


-ADq - Ds + ADt 


3730 ^"^^^^^^ 3^42^420 + ^75^220 
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